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The effects of weak point-like disorder on periodic systems at their upper critical dimension D c for 
disorder are studied. The systems studied range from simple elastic systems with D c = 4 to systems 
with long range interactions with D c — 2 and systems with D c — 3 such as the vortex lattice with 
dispersive elastic constants. These problems are studied using the Gaussian Variational method 
and the Functional Renormalisation Group. In all the cases studied we find a typical ultra-slow 
log log (a;) growth of the asymptotic displacement correlation function, resulting in nearly perfect 
transiationai order. Consequences for the Bragg glass phase of vortex lattices are discussed. 

I. INTRODUCTION 

The study of the physics of elastic systems in the presence of disorder has long been a subject of great theoretical and 
experimental interest, with applications to systems as diverse as charge density wavesa, Wigner crystalsa, magnetic 
bubbles&Q, etc. . These studies were given a further impetus with the discovery of high-Tc superconductors which 
were found to .be of the Type II kind. The resolution of the various theoretical challenges posed by these disordered 
vortex latticesO is crucial for the technological applications of high-T c materials. Early studies by LarkinB of the 
elastic system with point-like disorder indicated that disorder is relevant and that perfect long range transiationai 
order is destroyed below four dimensions. Since for a periodic system thermal effects dominate over that of disorder in 
dimensions d < 2, d — 4 and d — 2 are the upper and lower critical dimension for the relevance of disorder. However, 
the approximate model of independent random forces acting on each vortex line used fin. Ref. was insufficient to 
obtain the asymptotic properties of the disorder dominated phase. Scaling arguments™] suggested, within a purely 
elastic description, a slower, logarithmic growth of displacements. However, general arguments also tended to prove 
that disorder always favors the presence of dislocationsmL9 which thus invalidates the elastic approximation and leads 
to a destruction of the transiationai order. Although the precise nature of the disordered phase was a subject of much 
debate there was a thus a consensus that transiationai order was absent and that disorder leads to a glassy statepwith 
diverging barriers and pinningl30. Various experimental features such as the existence of a first order transition! 12 !! 13 !, 
presence of large dislocation free regions in decoration experiments^, and termination of the melting lineE3t3 could 
not be explained by these theories. nr - 1 

A quantitative theory of the elastic vortex lattice with point disorder was recently developecOO using both 
Functional Renormalization Group (FRG) and Gaussian Variational methods (GVM). These methods were originally 
introduced to study random manifold systemsEfl In the case of periodic systems both methods were found to 
be in qualitative and quantitative agreement, the FRGJjieing precise in d = 4 — e dimensions and the variational 
method being applicable for 2 < d < 4. These studiealjO yielded a description of the vortex lattice at all scales 
and demonstrated that while disorder produces algejp*Aie [growth of displacements at short length scales, periodicity 



changes the growth to logarithmic at large scalesB'BllaO'Eil thereby resulting in a power law decay of transiationai 
order. However, due to the abovementioned energy arguments indicating&lLj that disorder always generate dislocations, 
it was important to determine whether these results, derived in the elastic approximation would apply to realistic 
systems. These energy arguments, reexamined in Ref. [n] turned out to be incorrect and instead it was found that 
dislocations are unfavourable for weak disorder in d — 3. This led to the prediction of a distinct thermodynamically 
stable dislocation free glass phase called the Br agg g lass because of its nearly perfect transiationai order and perfect 
topological order. Furthermore it was proposecOEa that the phase seen experimentally at low fields was the Bragg 
glass, which accounted naturally for the firsjUptder transition and the decoration experiments. These predictions were 
confirmed by numericaEJcJ and analyticalc3~Ej studies. 

These early studies were made using dispersionless elastic constants, an approximation indeed valid to study the 
large scale properties of the system. However, in systems such as high T c superconductors the penetration depths can 
be large compared to the lattice constant, resulting in the possibility of non-local elastic constants at intermediate 
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length scales. The effect of such non local elasticity was qualitatively studied in Ref. g,|7| and in more detail recently 
in Ref. [34|,[55] where the Larkin length and the pinning force were computed, using a renormalization group procedure. 

In view of the physical importance and the remarkable properties of the Bragg glass phase, it is interesting to 
compute other physical properties, such as the positional order for intermediate length scales where the dispersion of 
the elastic constants is important. Intemstingly, this dispersivity of the elastic constants pushes the d = 3 Abrikosov 
lattice to its upper critical dimensionoE3. We thus address in this paper the more general question of the effects of 
disorder right at the upper critical dimension D c (for a discussion on what happens at the lower critical dimension 
see Ref. |l7| , ^6| -^9|) . Other examples of such systems which mimic the behaviour at d = 4 are single component 
elastic systems in d = 2 such as charge density waves where long range 1 Jr Coulomb interactions increase the effective 
dimension of the system to d = 4 i.e., D c = 2. Other motivations for this problem are theoretical. First, the 
logarithmic growth of displacements found for d < 4 has a prefactor proportional to e = 4 — d. It is thus interesting 
to study how the displacements grow at the upper critical dimension. Second, it also allows us to check and compare 
the results of the GVM and the FRG when the disorder is marginal. 

The plan of the paper is as follows. Section || is devoted to the study of simple isotropic elastic systems at the 
upper critical dimension d = 4 using both the GVM and the FRG. We compute the correlation functions and the 



characteristic length of the system, and compare the two methods. In Section III we study physical systems with 



dispersive elastic constants and anisotropy, in particular the d = 3 vortex lattice in the non local elastic regime. The 



conclusions are presented in Section IV 



II. ELASTIC SYSTEMS AT THE UPPER CRITICAL DIMENSION D c = 4 

A. The model 

We consider the simplest model of an isotropic elastic system where deviations from the equilibrium positions, 
characterised by the variables u a , are described in a coarse-grained manner by the following continuum hamiltonian 

H = J d d x[~(du a ) 2 + V(x)p(x)} (1) 

Here d is the spatial dimension and the index a runs over the number of components of the variable u. Depending 
on the physical origin of u, the above hamiltonian can be used to describe various disordered elastic systems. For 
example, the field u could denote the displacements from the equilibrium position for vortex lattices or could be a phase 
variable in the case of charge density waves, c is the non-dispersive elastic modulus and V(x) is the random potential 
describing the point-like impurities which couples to the density p. The random potential is gaussia n cor related 
(V(x)V(x')) = WS(x — x'). We use the decomposition of the density in terms of harmonics of the latticalSO 

p(x)~p (l-d a u a +J2 eiK(x ~ u) ) ( 2 ) 

where K are the vectors of the reciprocal lattice and po is the average density. Using (||) and the replica trick to 
average over disorder we obtain the following replicated hamiltonian 

H cS = fd d xJ2^(d<) 2 - J2 ^cos(K-(u a (x)-u b (x))) (3) 

Here T is the temperature, and irrelevant terms above the lower critical dimension have been discarded. To obtain 
the effects of disorder on all length scales one needs to study the theory retaining all the cosine terms in the repli- 
cated Hamiltonian (^). However, since we are mainly interested in the asymptotic physics which is governed by the 
fundamental K vector Kq, we retain only the lowest harmonic K = Kq in (|). We now study the hamiltonian (||) 
using both the GVM and the FRG. 

B. Gaussian Variational Method 

To solve (§), we follow the lines of Ref. |l6|Jilj , and introduce the variational Hamiltonian 
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H « = - 2 j (^2.< G ^ a S (4) 



The indices a, b denote the replicas and a, (3 denote the components of the field. Here G is the variational Green's 
function matrix which is parametrised as 

G- b ] a0 = f^(q)S ab -at (5) 

fapiq) — cq 2 $a[3 and the a are called the self-energy parameters. Since the disorder induced interaction in H e g is 
essentially local, the cr's are all taken to be constants. These parameters are determined by saddle point equations 
which are obtained by minimising the variational free energy F vaT = Fq + (-ff c ff — Hq)h with respect to a. The 
variational equations are (see Ref. |l7| for technical details) 

<t"(v) = K Koa K o0 e^ K ^ K ^t^ (6) 



where G^Jq) = E b G-l ab (q) and B is defined as 



There are two kinds of so lutio ns to these equations. One is the replica symmetric and the other is the replica symmetry 
broken one. Earlier workll3il3 has shown that the appropriate solution has a replica symmetry bxoken structure (RSB). 
We thus parametrize all quantities off-diagonal in the replica indices by hierarchical matriceaij 

G aa (q) = G(q) 

B aa {x) = B{x) (8) 
G ab (q) = G(q,v) 

where < v < 1 and B ab (x) = B(x,v). The indices a, (3 have been suppressed in the parametrisations presented 
above. 

The saddle point equations (||) now take the form 

a a0 (v) = K KoaKo0e l-hK Oa K o0 B^o, v)] (Q) 

where 

B af3 (0,v) = 2Tj -0- d [G aP {q) - G af3 {q,v)] (10) 

Since B(0 7 v) in (^) is essentially local, it follows that it is isotropic which results in the following simplification i.e., 
B a f}(x = 0) = 5 a pB(v) with B(v) defined below. Using this in (g) we find that cr Q ^(w) = 5 al 3a{v). We now search for 
a solution for a which has full replica symmetry breaking imbedded in it. This can be done in an easier manner by 
recasting the equations in terms of a new variable [a] = va(v) — J dua(u). .P.liigging in the ansatz of RSB, it can be 
shown that the correct solution for [a] has the following form for 2 < d < 4HEZL 

2 

[a] (v) oc vv , v < v c 

M(«)=E, v>v c (11) 

At the upper critical dimension there are corrections to the simple power law. Note that a knowledge of the exponent 
9, the break-point v c and £ is sufficient to fix the full functional form of tjae original cr(v). As a first step in evaluating 
these quantities one uses the the inversion rules for hierarchical matriceala to obtain 



B(0,v) =B(0,v c ) 



, , d d q 2Ta'(u) 

du / - — h 1 — (12) 
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where 

B(Q,v c )= [ ^- (13) 

Therefore, using Eqs.( pj],p^,^3[ ) in (^) one can obtain a self-consistent solution for a. Adapting the general framework 
presented above to the case of a single component model where Gj 1 — cq 2 and using ( |l2[ ) we find that the equation 
determining [a] is 

, , f d d q KqT 

^iw^-^MW]^ 1 (M) 



Solving for [a], we find 



[tr](v) = ^ A (15) 
1 log 2 ^ 



where A is the ultra-violet momentum cut-off and_A-j= K^TSd/Ici where 5^ is the value of the angular integration in 
d dimensions. For d < 4, [a] is a pure power lawEiltZl with the exponent 8 = d — 2 as shown in (jll]). The logarithmic 
corrections are specific to d = 4 (or rather d = D c as will be seen later) and have important consequences for the 
physical properties such as the correlation functions. The solution (|l5| ) for [a] is valid upto some v = v c above which 
[a] is a constant. This constant value is again fixed by (||). 
Using (Q) and (^) we can now calculate the correlation functions. The displacement correlation function is given by 



= 2T 
with 



B{x) = ((u(x) - u(0)) 2 ) (16) 
^(l-cos(g*))G (17) 



G {q ) = G c{q )[l+f% My,, J (18) 
Jo v [G c + [o](v)] 

Since the long- wavelength behaviour of B is governed by the small v behaviour of [a], the asymptotic displacement 
correlation function is 

B (x) = JLloglog(Az) (19) 



This growth of the displacements is indeed as expected slower than the logarithmic growth seen below four dimensions 
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B(x) = -^A d \og\x\ (20) 



with Ad = 4 — d and rn-is the number of components of u. 

One can also extractO the Larkin length R c given by R c — from the variational solution presented above. At 
zero temperature we obtain 



a 2c 2 tt 2 
2^ 6XP KqW 



^ = - ex Pl74TI7 (21) 



where a is the lattice spacing. The exponential growth, is characteristic of the upper critical dimension and is an 
extrapolation of the power law form 

R c = (cZ/WKfaY^-V (22) 

occurring for d < 4. The result ( |2l| ) given by the variational methnri-coincides with the naive dimensiDaal result or 
the more accurate result for R c obtained using the FRG procedureEjH Thus, as was already the caseESEJ for d < 4, 
the variational method gives the correct disorder dependence of the Larkin length. 
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C. Functional Renormalisation Group 



We now compute the correlation functions-, psing a functional renormalization group procedure. The results of the 
GVM and the FRG agree very well for d < 4l3'tZL As is usual with renormalization procedures, we can expect the FRG 
to work extremely well at the upper critical dimension. This provides another test to verify whether the equivalence 
between the FRG and GVM results which exists in d < 4 persists in d — 4. For simplicity, here we restrict to a single 
component model. 

We denote by A(u a (x) — u b (x)) the nonlinear term, proportional to the disorder strength in the hamiltonian (|^). 
For clarity we set Kq = 2tt. The function A(z) is initially ~ cos(27rz) and the idea is to obtain flow equations for 
this function whose fixed point solutions will allow one obtain the correlation functions. In dimensions d > 4, the 
gaussian fixed point A(z) = is perturbatively stable, and thus the cosine potential A can be expanded around z = 0. 
However, for d < 4, we need to retain all the terms in the expansion because they are all marginal. Hence, to study 
the flow of all these terms one uses the FRG as opposed to the usual RG. To obtain the flow equations of A and T, 
we use the rescaling x — * e l x and u — * e^u. Sipce the potential A is periodic in its argument, u cannot be rescaled, 
implying f = 0. Using a momentum-shell RGcS, we find that the system is described by a zero temperature fixed 
point. The corresponding flow equations for A and T for d < 4 is 

dT , 

^ = eA + i(A") 2 - A"A"(0) (23) 

where e = 4 — d and corrections to these equations are of higher order in e. A factor l/Sd = 2 d ~ 1 it d / 2 T[d/2] has been 
absorbed in A. From the flow equation for T we see that temperature is an irrelevant varia ble f or all d > 2. The 
periodic fixed point solution to this highly non-linear differential equation for A, for finite e islISIIZI 

A* = — (— — z 2 (l - zf) (24) 
72 V 36 V ' ' V ' 

for z in the primary interval [0, 1] and is extendable to other values of z by periodicity. The solution in ( p4| ) is 
non-analytic at z — and this non-analyticity is generic to disordered fixed points of the kind studied here. 

As anticipated, ( p4| ) shows that at the upper critical dimension the fixed point is A(z) = 0. Note that the RG 
equation in this case (d — 4) gives the exact large scale behaviour while in d = 3 one must rely on the O(e) expansion. 
To find the scale dependence of the function A;(z) in d = 4, we use the rescaling 

Af(z) = | (25) 

Substituting this in (E3l), the corresponding fixed point equation satisfied by y is 



o = y+(\y") 2 -y"y"(o) (26) 

Since this equation has a structure identical to ( |23| ) with e = 1 we can use the same solution to obtain 

»-s<5-'< 1 -""> (27) 



in the interval [0, 1]. Therefore, ( |27| ) together with the rescaling (|25|) for A* specifies the fixed point solution in d = 4. 
This solution also exhibits the afore-mentioned non-analytic behaviour at the origin. Another interesting feature of 
the solution to ( p3| ) in d = 4 is that the explicit I dependence of the solution shows us how the fixed point solution 
is approached as one scales the system. Using the above results, we can now calculate the displacement correlation 
function B(q) = (u a (q)u a (—q)} which satisfies the following flow equation 

B(q) = e il B(qe l ,T, A/) = {-f B{a~ l ,T = 0, A*) (28) 

qa 

In (|2§|), we have used e l — {qa)^ 1 where a is the lattice spacing which provides the ultra-violet cut-off. For long 
wavelength correlations, one can expand (pq) perturbatively in A* to finally obtain 

B(q)= (-)\-A*"(0)) (29) 



5 



Using ( p7| ) and restoring the factor S4 = 87r 2 this leads to 



1 

36 q A log(ag) 



B (l) = — 03 (30) 



The fourier transform of ( |3C| ) results in the following real space displacement correlation function 

8tt 2 



B(a;) = - u(0)) 2 ) = ^ log log Ax (31) 



where we have restored Kq. We find that both the GVM and FRG methods (|1_9|,^_1|) yield the same form for 

B{x) ~ -^logloglAx] (32) 

The prefactor b which also determines the exponent in the translation correlation function is b var — 2 and bf rg = 2.19. 
This discrepancy in the value of b is the same as the one that exists at dimensions d < 4. There such a difference was 
attributed to the fact that the variational method underestimates the effect of fluctuations and one expects the the 
same argument to hold in d = 4 as well. 



III. PHYSICAL SYSTEM WITH NON-LOCAL ELASTIC CONSTANTS 



A. Single Component Systems 

Let us now turn to a realistic system with long range interactions. The simplest system is a CDW system with 
long range Coulomb interactions. In this case the field u is a simple scalax-field describing the phase. Due to the long 
range nature of the Coulomb interaction the elastic propagator becomesEll 

G- 1 {q)^cq 2 +dq (33) 

instead of the simple q 2 propagator corresponding to local elasticity. By writing dq = a(q)q 2 with a = |, we can view 
the Coulomb contribution as being equivalent to a non-local elastic modulus. When d is large this shifts the "physical" 
upper critical dimension from D c — 4 to D c = 2. The results of the previous section are thus directly applicable. We 
find that [a] has the same form as that given by @ with A= 2 f . Rather surprisingly, the asymptotic displacement 
correlation function is found to grow in precisely the same manner as the correlations in d = 4 

B( 2 ;) = Aioglog(A a ;) (34) 

The FRG can be used in the present case with the replacements e = 2 — d and £4 — > S2 in |3(]). This results in 
the same correlation as in ([34]) with the prefactor 2.19 instead of 2. Note however, that in systems such as Wigner 
crystals where transverse displacements exist, these transverse displacements-ace insensitive to the long range part of 
the Coulomb interaction and the upper critical dimension remains D c — 4 c!Jcd, leading to a different result for the 
correlations in d = 2. 



B. Vortex Lattice 



We focus in this section on another physical realization of elastic systems at the upper critical dimension, namely, 
the vortex lattice in d = 3. The physical properties of such a system with non-dispersive elastic constants are now 
well known (see e.g. Ref. [37] and references therein) and lead to the phase diagram of Figure 
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FIG. 1. Schematic phase diagram of the 3-d vortex lattice in the H — T planecJ'L-J. The stability region of the Bragg glass 
phase in the magnetic field H , temperature T plane is shown for moderate fields (the very low field region is not shown). The 
thick line is expected to be first order, whereas the dotted line should be either second order or a crossover. Upon increasing 
disorder the field induced melting occurs for lower fields as indicated by the thin solid line. 



The Bragg glass occurs for weak disorder where collective pinning works. The system is dislocation free and 
displacements are logarithmic. At high fields the disorder is stronger and dislocations are present. We study here the 
effect of non-local elastic constants in the Bragg glass phase. For a realistic form of the dispersion of these moduli i.e., 
momentum dependence, the system could effectively behave like a d = 4 system at intermediate lengthscales. This 
scenario is realised in high-T c vortex lattices. 



C. Model 



The prototype elastic hamiltonian for d = 3 vortex lattices formed in type II superconductors with point-like 
disorder is 



H 



d 2 q±dq 2 



[(en - c 6e )(q a u a ) 2 + c 66 {q a upY 



(2tt) 3 

+ c i4: (q z u a ) 2 ] + / d 2 rdzV(r, z)p(r, z) 



(35) 



Here u are the displacement vectors and a, (3 represent the x,y components. In (|3q), Cn, Cqq and C44 are the 
compression, shear and tilt moduli respectively. These elastic constants are in general momentum dependent and 
typically depend on the detailed properties of the superconductors. They have been evaluated within a Ginzburg 
Landau theory and are given by the following expressions^ 



C66 



C44 



Cll 



(8ir\) 2 
B 2 



1 



4tt 1 + \ 2 /e 2 q 1 2 + \ 2 q 2 z 
B 2 



-44 



(36) 



1 + \ 2 /e 2 q 2 



4tt (1 + AV)(1 + A 2 /eV 2 + A 2 g2) 



Here $0 is the flux quantum, B is the magnetic field along the c axis, e = J M± /M z C 1 is the anisotropy 
parameter, A is the London penetration depth in the ab plane and q± is the in-plane momentum. c| 4 is the single 
vortex contribution which is negligible for q± < 1/a. Note that though the shear modulus is always non-dispersive, cn 
and C44 are both dispersive, in other words, non-local. The origiMj-oLt liis dispersion is known to lie in the long ranged 
Yukawa like interaction between vortex segments. Earlier workBQ'taO'cS had focussed on regimes of the field where 
the dispersion in these constants was negligible and one recovered a completely local hamiltonian. A logarithmic 
growth of the asymptotic transverse and longitudinal displacements was found in this case. This growth was also 
found to be isotropic despite the obvious anisotropy in the elastic hamiltonian. Apart from this logarithmic regime, a 
random manifold regime with a power law growth of displacements for intermediate length scales where all the cosine 
terms were relevant was also found. However, in addition to these regimes, there exists another regime where the 
dispersive nature of the elastic constants discussed above becomes crucial. This is the region where cn and Cqq are 
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both non-dispersive with cu» cqq and C44 is dispersive (or en is dispersive but much larger than cqq in which case 
we can set c^ 1 = in the following). The approximation to C44 in this regime is given by 



B 2 e 2 

C44 = i^qj (37) 

This q— dependent elastic modulus changes the nature of the effective long wave-length theory. In Ref. 16,17], it 
was argued that such a non-local regime might manifest itself in cases where the translational correlati on l ength R a 
was smaller than j. The effect of such non-local elasticity on vortex lattices was reexamined recentlyaEJ, and the 
disorder dependence of the Larkin length was computed using a second order FRG calculation. We focus here on the 
calculation of the translational order parameter. As will be shown below, correlation functions in this regime behave 
very differently from that of the logarithmic and random manifold regimes. 

The disorder independent part of the Hamiltonian when rewritten in terms of the transverse and longitudinal 
components in Fourier space is 

H = [ - ,V~ff Z (C66ffl + C44i?z)«T + (en?! + c 44<?f )u\ (38) 



., (2tt) 3 

Using the formalism presented in Sec. II and the following connected Green's functions 

= cuql + c 66 q 2 ± 

G~l = c m 2 z + Cll ql (39) 

we find that [a] in the present case is determined by the following transcendental equation (refer to the Appendix for 
details) 



■ n 2 A 2 r -3/2 . -3/2, 647rVf 

M lQ g p[ C H + C 66 ] = K 2 T V ( 40 ) 

where T — B 2 e 2 / (inX 2 ) and A = — is the lattice cut-off. Solving (J4C|) , we find that [a] has the same form as that 
given in ([15]), with the constant A now given by 

647T /=. 1 1 . i . 

^ = + — ] (41) 

K ° T C?! C| 6 

Using the result for [a] the displacement correlation functions are given by 

l 2 q_i_a 



B L , T (r,0) = 2T I [cos 2 (g)(l - cos(q ± r co S (6)))G L<T (42) 



+ (1 - cos 2 (60)(l - cos(<ur cos(6»)))G t ,l] (43) 



Here 



/=. / \ f 1 dv MM „ .... 
G L , T (q ± ,q z )= - G cL , T (44) 

Jo v G cL T + [(t\ 

Performing a whole lot of cumbersome integrations (see Appendix), we finally obtain the following leading terms 
in the displacement correlation functions: 

B L>R (r, 0) = B(r, 0) = log log A V (45) 

Note that this behaviour of the correlations is exactly analogous to the ones seen earlier in d = 4 (cf. jl9|) , (|3l|) ) and 
d = 2 proving again that non-local elastic constants of the type considered here increase the effective dimensionality of 
the system. Note that despite the strong anisotropy in the hamiltonian, isotropy is restored in the the long wavelength 
correlation functions in this dispersive regime too. Another interesting point is that these correlations are universal in 
that they do not depend on the value of any parameter in the theory. The above result can now be used to calculate 
the translational order parameter. Within the Gaussian approximation it is thus given by: 
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T = (cxpiK Q .(u(r)-u(0))) = exp-^((u(r, 0) - u(0, O)) 2 ) (46) 

* (m7) ^ 

Note that the decay of Or in the present case is much slower than the - decay seen in the case where the elastic 
constants are local. This implies that the system will look more ordered than it does in the other regimes. 



D. Physical Discussion 



We now estimate the regime of validity of the previous calculation i.e., the range in which the dispersive approxi- 
mation to C44 and hence the log log growth of displacements are valid (for a very complete description of the possible 
regimes in anisotropic superconductors see Ref. |35j). Re-analysing the expressions ( pq) for the elastic constants we 
see that the approximation ([37]) to C44 is valid provided one is working at length scales L± and L z which satisfy the 
following bounds 

L ± < A/e (48) 
L x « L z /e (49) 

It is also interesting to check when the above dispersive regime includes both the Larkin lengths R± and R z . The 
Larkin lengths are determined by the balance between the elastic energy and the disorder which in d = 3 takes the 
form 

c 66 r 2 f R z 2 ~ c 4i rjR- 2 ~ r f R ± Rl /2 W 1/2 (50) 



Here rj is the length over which disorder is correlated. We immediately see that R z /R±_ = \J c^/cqq and using (|3 
we thus obtain 

R z = eR\ja (51) 

This in conjunction with ( f48|) yields the following bounds on R± 

a < R ± < A/e (52) 

The upper limit corresponds to R z — X 2 /(ea) and the lower limit R z = ea. The latter may be smaller than the 
distance d between layers resulting in complicated dimensional crossovers. To stay away from possible dimensional 
crossover we will restrict R± such that 



V^d/e < R± < A/e (53) 

which corresponds to d < R z < A 2 /(ea). Note that the above regime becomes non-existent above the "decoupling 
crossover field" B cr — <& e 2 /d 2 . 



B(r) t 










Larkin 


Dispersive 
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log(r/a) 

FIG. 2. The growth of the displacements B(r) as a function of log(r/a) in various regimes. 
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We therefore see that the result of ([l5|) can be used at all lengths scales within the bounds prescribed by ([l8|) . For 
asymptotically large scales, the elastic constants become non-dispersive and one expects a crossover to the simple 
logarithmic growth of displacements discussed in Ref. Il6],[l7|. These results are summarised in Fig. 2 . 



IV. CONCLUSIONS 



We have studied the long wavelength correlations of elastic systems at their upper critical dimension in the presence 
of weak pointlike disorder. The cases analysed include the simple elastic system at d = 4, single component systems 
with i interactions and vortex lattices with non-local elastic constants. The systems were studied using the GVM 
and the equivalence between the GVM and the FRG methods was explicitly demonstrated in the simpler cases of 
d = 4 and d = 2 and we expect such a correspondence to hold in the other cases as well. In all the cases studied we 
observe the striking fact that the growth of asymptotic displacements is generic in that it is always loglog(a:). This 
was found to be true even in the case of the anisotropic vortex lattice, where isotropy was recovered in the asymptotic 
correlation functions. In addition, the prefactor of this log log term (which determines the decay of the translational 
order within the GVM) is independent of the parameters of the models considered. These systems are again Bragg 
glasses and the peaks in the diffraction will be more pronounced because the systems are more ordered than their 
other disordered counterparts which typically show a logarithmic growth of displacements. This feature should be 
verifiable using decoration and diffraction experiments. 



V. APPENDIX 



We first write down the explicit form of the equation which determines [a] in the case of the vortex lattice described 
in Sec. III. This is a generalisation of ( |l4| ) to the two component case and it reads 

.KpT r d d q l l 

a[V) 2 J (2nr [ (G cT (q)- 1 + W](v)) 2 + (G cL (g)-i + H(«))2 J - {b > 



Using the expressions for G c t,l given in (|39|) and integrating over q± and q z we arrive at (|40|). 

Next we present the details of the calculation of the displacement correlation function discussed in Sec.V . Performing 
the angular integrations in Qwe get 

Bl = / q±dq±dq z [(l - J (q±r) + ^^)G L (55) 

(2-Ky J 2 q±r 

1 Ji{q±_r) - 

+ )\G T (56) 

2 q±_r 

Jo and Ji are the usual Bessel functions and the expressions for Gl,t are given in (^i|). ( [561) is rewritten as 

2T r dv 
(2^J V 2 



Bl = —\ ^(h+h) (57) 



where 



and 



h= f q±dq ± dq z [ ^ 4 ]f „ V [l2 ](^ - J° + ~) (58) 
J Fit + CnQl + cnq\ + [a\q± 2 q±_r 



l Tqp + ceeq^Tq? + c m q\ + [a]q 2 ± iK 2 q ± V 



h = I q±dq ± dq z [ ^ 2n ]\ 2 ^ - ^ - . . 2 ](- - ^) (59) 



Integrating over q z we get 



h = ^= /tAt ^l!! 1 "^ 1 ]^ - Jo + A) (60 ) 
2VT^i y± qM y/ Cll qi + [a] U 2 q^r* K ' 
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li has the same form with c\\ «-> cqq and the Bessel functions appropriately changed. Integrating with respect to q± 
we now get 



2r y en 

" -(/ (a)Ko(a)-/i(a)^i(a))] 



4ci 



(61) 

(62) 
(63) 



and 



h = — ; o + og 2A, 

2y/T^ 1 r 2 4c 66 S V M 



1 ff 



2r V c 66 



(Jo(/3)^i(/3)-/i(/3)if (/3)) 



(64) 
(65) 



Here A is the momentum cut-off and the / and K are the modified Bessel functions whose arguments are a = £\/ — 
and P = 7j\ — • Finally, adding up the contributions /i and l<i we get 



Si(r,0) 



47rVrJ « 



dv 



(66) 



The transverse correlation function Bt is obtained from the above expression by interchanging en and cqq. Using 
the solution for [a] given by (Eoh and 



Br 



A V T f dv A v v 2 E r 2 A 



-[lo. 



log- 



327rVT 7 w 2 L b A 2 J 1 2 * A v v 



CE 



(67) 



3/2 3/2 

where E — c xx -h c 66 . Performing the f integration gives us the result 



B i)T (r,0)~-jloglogAV 



(68) 
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